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ABSTRACT 



43 

Context. Recent spectropolarimetric observations of Herbig Ae/Be stars have yielded new arguments in favour of a fossil origin for 
j_( . the magnetic fields of intermediate mass stars. 

Aims. To study the evolution of these magnetic fields, and their impact on the evolution of the angular momentum of these stars during 
' the pre-main sequence phase, we observed Herbig Ae/Be members of young open clusters of various ages. 

Methods. We obtained high-resolution spectropolarimetric observations of Herbig Ae/Be stars belonging to the young open clusters 
NGC 661 1 (< 6 Myr), NGC 2244 (~ 1.9 Myr), and NGC 2264 (~ 8 Myr), using ESPaDOnS at the Canada-France-Hawaii Telescope. 
Results. Here we report the discovery of strong magnetic fields in two massive pre-main sequence cluster stars. We detected, for the 
first time, a magnetic field in a pre-main sequence rapid rotator: the 10.2 M Q Herbig B1.5e star W601 (NGC 6611; v sin i ^ 190 km/s). 
Our spectropolarimetric observations yield a longitudinal magnetic field larger than 1 kG, and imply a rotational period shorter than 
1.7 days. The spectrum of this very young object (age ~ 0.017 Myr) shows strong and variable lines of He and Si. We also detected 
a magnetic field in the 12.1 M e Bl star OI 201 (NGC 2244; v sin i - 23.5 km/s). The Stokes V profile of this star does not vary over 
5 days, suggesting a long rotational period, a pole-on orientation, or aligned magnetic and rotation axes. OI 201 is situtated near the 
, Zero-Age Main Sequence on the HR diagram, and exhibits normal chemical abundances and no spectrum variability. 
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1 . Introduction Their rotational periods are typically about 3 days, whereas nor- 
mal A/B stars rotate with periods less that one day. Remarkably, 

About 5% of main sequence (MS) A and B stars host strong, some Ap/Bp starfj haye significantly i onger rotational periods, 

organised magnetic fields. All these magnetic stars also show which in extreme cases can approach one century . Stepien 

*3 . peculiar and variable chemical abundances, and are called the (2000) considered different hypothesis to understand the origin 

magnetic Ap/Bp stars. The origin of these magnetic fields is of the slow rotation of Ap/Bp stars . He concluded that only mag- 

still not well understood. The favoured hypothesis is the fossil netic brakingi acting during the PM S phase of stellar evolution, 

field model, which proposes that the magnetic fields are relics couM explain the slow rotation of the Ap/Bp stars Tf Stepie n's 

of fields present in the molecular clouds from which these stars conc i usion is correct then some HAeBe stars should host both 

formed. This implies that some pre-main sequence (PMS) stars magnetic fields and disks, and an evolution of the rotation of 

of intermediate mass should also be magnetic, with field prop- thefje stars shouW occur duri their pMS phase 

erties similar to those of the Ap/Bp stars. Many such PMS stars „ ,. „ „ £T ta ™ 

, . „ „ i • a m : ti n , Our recent studies of the magnetic fields of HAeBe stars 

are identifiable observationally as Herbig Ae/Be (HAeBe) stars. , , , . ,. c , .. n . „ , , . 

J ' have led to the discovery of 4 new magnetic HAeBe stars, bring- 
Furthermore, most of the Ap/Bp stars rotate very slowly j stro ng, new argum e nts in favour of th e fossil-fie ld hypoth- 
compared to normal (nonmagnetic and nonpeculiar) A/B stars. gsis d Wa deetal.l S ICatala et al.1 120071: lAlecian et all [2008. 
Folsom et al. in prep.). 

Send offprint requests to: E. Alecian — jr^u u i_ ,,j»u i t c 

1 To proceed further, we have begun to study the evolution of 

Based on observations obtained at the Canada-France-Hawaii .• j..- c ■ * j- . u u 

_ , ._„.„_. .... , . . „ , „ .. magnetism and rotation of intermediate-mass stars by observ- 

Telescope (CFHT) which is operated by the National Research Council TT . _ . , n ,,. 

of Canada, the Institut National des Sciences de 1' Univers of the Centre ln 8 HAeBe stars ln y oun § °P en clusters of dlfferent a § es - We 

National de la Recherche Scientifique of France, and the University of have obtained ESPaDOnS spectropolarimetric observations of a 

H awa ii sample of 61 young intermediate-mass stars in the open clusters 

** visiting astronomer at the Pic du Midi Telescope Bernard Lyot NGC 2244, NGC 2264, and NGC 661 1. While a full description 

(USR5026) of Universite Paul Sabatier and CNRS Institut National des of the experiment and the results will be provided in a forth- 

Sciences de 1' Univers coming paper, in this Letter we report the detection of magnetic 
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Table 1: Log of observations of W601 and OI 201. Columns 1 
and 2 give the UT date and the Heliocentric Julian Date of the 
observations. Column 3 gives the total exposure time. Column 
4 gives the peak signal to noise ratio per 1.8 km/s CCD pixel 
(at ~ 708 nm for W601, and at ~ 552 nm for OI 201) in the 
spectra, and Column 5 gives the signal to noise ratio per 1.8 km/s 
in the LSD Stokes V profiles. Column 6 gives the longitudinal 
magnetic field. 



Date (d/m/y) 


HJD 




S/N 


S/N 


Be 


UT Time 


(2 450 000+) 


(s) 




(LSD) 


(G) 




W601 








10/08/06 6:09 


3957.76004 


3600 


220 


1910 


-514+335 


06/03/07 14:28 


4166.10187 


3600 


180 


1610 


-313+427 


06/03/07 15:35 


4166.14879 


3600 


120 


1010 


-1279+696 


09/03/07 14:34 


4169.10634 


3600 


190 


1760 


1391+394 




OI 201 








04/03/07 9:49 


4163.91136 


960 


150 


1600 


-477+73 


05/03/07 6:12 


4164.76045 


960 


160 


1820 


-543+65 


06/03/07 6:47 


4165.78497 


960 


160 


1710 


-546+68 


09/03/07 9:03 


4168.87910 


960 


150 


1490 


-471+76 



fields in two young, massive stars: W601 (= NGC 661 1 601) and 
OI 201 (=NGC 2244 201). 

2. Observations 

We observed W601 and OI 201 in August 2006 and March 2007, 
using the ESPaDOnS spectropolarimeter installed at the Canada- 
France-Hawaii Telescope (CFHT, Hawaii). The observing and 
reduction procedur es are fully described in many papers (e.g., 
Alecian et al. 2008) and provide us with the intensity / and cir- 
cular polarisation V spectra of a star, with a resolving power of 
about 65000 and S/N ratios ranging from 120 to 220 per CCD 
pixel for W601, and from 140 to 160 per CCD pixel for OI 201. 
The Libre-Esprit reduction package, developed especially for 
ESPaDOnS, calculates, i n addition to the I a nd V spectra, a null 
polarisation spectrum N dDonati et al.ll 1997b . which allows us to 
check that the signature in the Stokes V spectrum is not spurious. 
The log of the observations is reported in Table 1 . 

To increase the effective signal to noise (S/N) ratio of 
our data, we appl ied the Least Squares Deconvolution (LSD; 
iDonati et al.ll 1997b procedure using tailored line masks of 500 
and 680 lines for W601 and OI 201, respectively. The S/N ratios 
of the LSD Stokes V profiles are at least 10 times larger than the 
S/N ratio in the original spectra (Table[TJ, and allow us to clearly 
see strong Zeeman signatures in each of our observations. 

3. W601 (NGC 6611) 

The B1.5 star W601 is member of the open cluste r NGC 661 1, 
situat ed at a distance of 1.8 kpc from the su n dDufton et al. 
2006), and has been classified as a PMS star by iMartavan et al. 
(2007). Due to scatter in the colour-mag nitude diagram, t he age 
of the cluster is not well constrained. Ide Winter et al.l (1997) 
conclude that the sample of PMS stars is not coeval, with ages 
spread from to 6 Myr. According to the coordinates of the cen- 
ter of the cluster of iBelikov et all(11999l) the position of W601 
is offset, while still situated inside the limi ts of the cluster . Its 
proper motion from the Tycho- 2 catalogue (H0g et al. 200CJ), its 
reddening law (Guarcello et al. 2007), and its position in the HR 
diagram dDufton et al.l2006l) are fully consistent with NGC 661 1 
membership. 




4.45 4.40 4.35 4.30 4.25 

log(Teff) 



Fig. 1: A pre-main sequence Hertzsprung-Russell diagram 
showing the positions of W601 (red square) and OI 201 (green 
diamond). The CESAM PMS evolutionary tracks (black full 
line), the zero-age main sequence (ZAMS) (dot-dashed line), 
and the birthline (blue dashed line) are also shown. 



According to [Kumar et al.l (120041) . this star exhibits mid- 
infrared (MIR) excess. In our data, we observe emission in Ha 
and H/3. Given its membership in NGC 66 1 1 , its HR diagram 
position, and its spectral properties, we concur with the classifi- 
cation of W601 as a PMS Herbig Be star. 

Using the effective temperature (lo g T e g = 4.35 + 0.05) and 
luminosity (log L/L Q = 3.96 + 0.11) of D ufton et ail d2006l) . we 
compared the HR diagram position of W601 with PMS evolu- 
tionar y tracks of so lar abundances, calculated using the CESAM 
code dMorelll 1997th We find a mass M = 10.2+'^ M G and a ra- 
dius R = 6At.fl R e>' leadin g to lo g,? = 3 - 8 -o3 ( C S S )- The a g e 
of the star has been estimated to be about 0.016^ ^3 Myr from 
the birthline ca lculated with a proto stellar mass accretion rate of 
10- 4 Mo.yr- 1 dPalla & Stahlerl 19921) (Fig.Q]). 

We obtained spectra of W601 sampling a variety of 
timescales, from one hour to several months (Table [TJ. We ob- 
serve emission in the cores of the Ha and H/3 Balmer lines, and 
in the wings of Ha, varying on a timescale of one hour. The 
absorption lines are strongly broadened, with vsin/ = 190 ± 
10 km.s -1 . According to our determination of the stellar radius, 
this implies that the rotation period of W601 is shorter than 1.7 d. 

Most of the spectrum of W601 is consistent with a syn- 
thetic spectrum of T e ff = 22500 K, logg = 4 (cgs), broad- 
ened by vsin/ = 190 km.s -1 , calculated using TLUSTY non- 
LTE atmosphere mode ls and the SYNSPEC code dHubenvl 19881: 
iHubenv & Lanz|[l992l) . However, we observe that all He 1 lines 
are substantially stronger than the synthetic ones calculated us- 
ing the solar abundance (Fig. [2] lower frame). The He lines 
alos show variability in their depth and shape on a timescale 
of 1 hour. We also observe strong and variable Si m lines, 
but with a shape variability that differs from that of the He 1 
lines. These characteristics suggest that W601 is a PMS He- 
strong star with He and Si ab undance spots on its surface (e.g., 
Bohlender & Lands treelll990l) . The observed variations of the 
He and Si lines on a timescale of one hour (which we interpret 
as rotational modulation by the abundance spots) are consistent 
with our predictions of a rotation period shorter than 1 .7 d. The 
line profiles show no systematic, radial velocity variability sug- 
gestive of binarity. All the LSD / profiles and Stokes / spectra 
obtained on the 3 nights show the same small structures in the 
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Fig. 2: Upper frame: LSD / (bottom), V (upper) and N (mid- 
dle) profiles of W601. Full black line: 10/08, red-dashed line: 
6/03 14:28, green dot-dot-dot-dashed line: 9/03. The mean error 
bars are plotted next to each profile. Lower frame: The broad, 
variable He i /15015 of W601, observed on 6/03 (full-black line) 
and 9/03 (red-dashed line). The green long-dashed-line is a 
TLUSTY/SYNSPEC synthetic profile computed for a solar He 
abundance. 



line profiles. Their origin is unknown, although they are appar- 
ently not instrumental. 

We plot the LSD Stokes / and V profiles obtained on three 
different dates in Fig. [2] (upper frame). We observe some vari- 
ations in the shape of the / profiles. The red side of the profile 
has changed between August 2006 and March 2007, while the 
blue side varied slightly between March 5th and March 8th. The 
mask used to calculate these profiles contains mainly He i lines. 
Therefore, the variations in the Stokes / profile reflect primarily 
those observed in the He i lines of the Stokes / spectrum (Fig. 
[2] lower frame). The polarisation Stokes V profiles, plotted at 
the top of Fig |2l show strong and variable Zeeman signatures as 
broad as the photospheric / profile, and centred on its centroid, 
which indicates that this signature is created by a magnetic field 
at the surface of the star. The comparison of the LSD V (upper) 
to the LSD N (middle) profiles supports that the V signature is 
real, and not a spurious signature produced by the instrument or 
data reduction procedure. 

The LSD V profiles also vary on a timescale of one hour. 
Our previous studies of magnetic HAeBe stars have shown that, 
like the Ap/Bp stars, they host large-scale organised magnetic 
fields and satisfy the oblique rotator model (Alecian et al. 2008; 



Folsom etal., in prep.). This model basically consists of a dipole 
placed inside the star , with a magnetic axis inclined with respect 
to the rotation axis dLandstreetl 1 1 970h . As the star rotates, the 
configuration of the surface magnetic field seen by an observer 
changes, and the observed polarisation signal therefore varies. 
Most Ap/Bp and magnetic HAeBe stars show temporal varia- 
tions of their polarisation spectra, which are modulated accord- 
ing to the rotation period of the star. By analogy, we propose 
that the variations of the Stokes V profile of W601 also reflect 
the rotation of the star. 

Us ing the Stokes / and V profiles and Eq. (1) of IWade et al.l 
(2000), we measured the longitudinal magnetic field B( of our 
data. The values of B( are listed in Table 1 and are found to 

vary between 1300 and ~ 1400 G, indicating that the mean 

magnetic field modulus is probably larger than 3 kG. 

Our data are insufficient to allow us to derive the oblique ro- 
tator model parameters of W601 . Additional spectropolarimetric 
observations with dense sampling of the rotation cycle are re- 
quired to accurately determine the rotation period as well as to 
model the intensity and the topology of the magnetic field. 

4. OI 201 (NGC 2244) 

The Bl star OI 201 belongs to the cluster NGC 2244, ilu- 
minating the Rosette Nebula. The cluster is located at a dis- 
tance of 1.7 kpc from the sun, with an age of 2.3 + 0.2 Myr 
dHensberge et al. 20001). The proper motion and radial veloc- 
ity of OI 201 (IChen et alJ l2007l: iPerezll 19911) . its reddening law 
dDucafi et alj|2003l) . its position in the HR diagram, and its ge- 
ographical position in the central blister of the cluster confirm 

unambiguously that it is a member of N GC 2244. 

OI 201 does not show an IR excess dBalog et al. 2007) and 



has not been detected as a n X-ra y source ( Wang et al. 2007b . 



According to iPark & Su ng ( 2002), this star does not have Ha 
in emission, which we confirm with our data. The position of 
OI 201 in the HR diagram is above, but close to the ZAMS - 
therefore that star is either an old PMS star or a young MS star. In 
both cases, this star is very young, but more evolved that W601. 

We used the photometric data of IPark & Sund d2002l) to de- 
termine the luminosity of the star: log(L/L Q ) = 4.05 + 0.08, 
while its temperature has been determined by IVrancken et al.l 
d!997l) : Jeff = 27300 ± 1000 K. We compared the HR dia- 
gram position of OI 201 with CESAM PMS evolutionary tracks. 
We find M = 12.1+^ M and R = 4.7^ M , leading to 
lo gg = 4.18 + 0.06 (cgs), consistent with the spectroscopic value 
of lVrancken et alj dl997). Assuming that OI 201 is PMS, its age 
since the birthline is estimated to be 0.041+°°Sg M Y r ( Fi §- IB, 
which is inconsistent with the age of the cluster. Therefore this 
star is more likely a MS star with an age of 3*2 Myr, estimated 
using CESAM MS evolutionary tracks. 

The four spectra, all obtained during a single run of 5 days, 
show no evidence of emission, either in the Balmer lines or 
in the He i and O i lines often observed in emission in PMS 
stars. We observe no variability in any lines and we find that a 
TLUSTY/SYNSPEC synthetic spectrum with r eff = 25000 K, 
logg = 4.25 (cgs) and solar abundances, broadened with 
v sin i = 23.5 + 0.5 km.s~', reproduces well the entire spectrum 
of OI 201 (Fig.[3]lower frame). Figure[3](upper frame) shows the 
LSD profiles of OI 201 observed on 4 different dates. The four 
Stokes / profiles are perfectly superimposed, illustrating the lack 
of variability in the spectrum. No structures are visible in the line 
profiles similar to those detected in W601. 

The Stokes V polarisation profile shows a strong Zeeman 
signature as broad as the / profile, and centred on its centroid, 
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Fig. 3: Upper frame: LSD / (bottom), V (upper) and N (middle) 
profiles of OI 201. Full black line: 4/03, red-dashed line: 5/03 
, blue dot-dashed line: 6/03, and green dot-dot-dot-dashed line: 
9/03. The mean error bars are plotted next to each profile. Lower 
frame: The sharp, nonvariable He i /15015 of OI 201, observed on 
4/03 (full-black line) and 9/03 (red-dashed line). The green long- 
dashed-line is a TLUSTY/S YNSPEC synthetic profile computed 
for a solar He abundance. 



which indicates that this signature is created by a magnetic field 
with an important poloidal component at the surface of the star. 
We note that the LSD N profile is totally flat compared to the 
Stokes V profile, indicating that the Stokes V signature is real. 

The Stokes V profiles of the four spectra observed between 
March 3rd and March 8th have not varied, and yield similar 
measurements of the longitudinal magnetic field found around 
-500 G (Table Q]). According to the oblique rotator model de- 
scribed in Sect. 3, the lack of variation in the Stokes V profiles 
indicates that either the rotation period of the star is very long, 
the rotation axis of the star is nearly parallel to our line-of-sight, 
or the magnetic axis is nearly aligned with the rotation axis. With 
our determinations of v sin i and R, we find that the maximum ro- 
tation period of OI 201 is 10.1 d. Therefore a very long period 
(-months) is not consistent with the vsin/ and radius. On the 
other hand, if the inclination of the rotation axis i to the line-of- 
sight is very small (say lower than 10°), rotational modulation 
would be very weak, and perhaps not detectable. Such a pole- 
on scenario would also imply a rotation period similar to that 
of W601 (for i = 10°, v eq = 135 km/s and P mt = 1.8 d). On 
the other hand, such a small value of i is statistically very un- 
likely (less than 1%). Furthermore OI 201 is the second star of 



our Herbig survey that shows no variability in the V profile (the 
first one being HD 190073: ICatala et al . 1 120071) . The probability 
to observe two pole-on stars in a same sample is small. A final 
option is that OI 201 has its magnetic axis aligned with the rota- 
tion a xis, which has been ob served in some MS He-strong stars 
(e.g jBohlender et al.ll 19871) . Additional spectropolarimetric ob- 
servations are required to distinguish between these different ge- 
ometries. 



5. Discussion 

We have discovered strong magnetic fields at the surfaces of two 
hot, very young massive stars, both of which are members of 
young open clusters. Although these stars have similar masses 
and are both magnetic, they show important physical differences: 
W601 is a rapid rotator, exhibits strong and variable lines of He 
and Si, and is situtated in the first half of the PMS phase. OI 201 
appears to be a slow rotator, shows normal photospheric chem- 
istry, and may already be on the MS. 

The presence of chemical peculiarities in the atmosphere of 
W601 demonstrates for the first time that such anomalies can 
develop very early during the PMS evolution of a B star. This 
prompts the question: why are no similar chemical anomalies 
observed in the atmosphere of OI 201? This star is only about 
20% hotter and more massive than W601, and it hosts a similar 
strong magnetic field. Moreover, its slower rotation and apparent 
lack of a circumstellar accretion disc (as suggested by its lack of 
emission lines and IR excess) would suggest that its atmosphere 
is likely more stable than that of W601. We suspect that the dif- 
ferences in the photospheric properties of these stars will yield 
important new constraints on models of chemical fractionation 
in the atmospheres and winds of B-type stars. 

Unless the rotation axis inclination of OI 201 is smaller than 
about 10°, it certainly rotates much more slowly than W601 and 
other B-type stars. Could its advanced age, interpreted within the 
context of magnetic braking, explain the large difference with 
respect to W601? 

The answers to these questions must await the acquisition 
of additional spectropolarimetric data. However, it is clear from 
the observations presented here that W601 and OI 201 repre- 
sent two remarkable B-type PMS stars that have the potential 
to contribute considerably to our understanding of magnetism, 
rotation, and evolution of intermediate-mass stars. 
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